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ABSTRACT 

Weak gravitational lensing surveys are rapidly becoming important tools to probe 
directly the mass density fluctuations in the universe and its background dynamics. 
Earlier studies have shown that it is possible to model the statistics of the convergence 
field on small angular scales by suitably modeling the statistics of the underlying den- 
sity field in the highly non-linear regime. We extend such methods to model the 
complete probability distribution function, PDF, of the shear as a function of smooth- 
ing angle. Our model relics on a simple hierarchical Ansatz for the behaviour of the 
higher-order correlations in the density field. We compare our predictions with the re- 
sults of numerical simulations and find excellent agreement for different cosmological 
scenarios. Our method provides a new way to study the evolution of non-Gaussianity 
in gravitational clustering and should help to break the degeneracies in parameter 
estimation based on analysis of the power spectrum alone. 

Key words: Cosmology: theory - gravitational lensing - large-scale structure of 
Universe Methods: analytical - Methods: statistical -Methods: numerical 
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1 INTRODUCTION 

Recent detection of weak lensing (see, e.g., Bartelmann & 
Schneider, 2001, for a review) due to the large-scale struc- 
ture in the universe has presented a new opportunity to 
study cosmology (see, e.g., Bacon, Refregier & Ellis, 2000, 
Hoekstra et al., 2002, Van Waerbeke et al., 2000, and Van 
Waerbeke et al., 2002). With a tight control on various sys- 
tematics these studies were able to relate the weak distor- 
tions in the images of high-redshift galaxies with the in- 
tervening large-scale mass density inhomogeneities in the 
universe. Such studies which are statistical in nature (as 
opposed to strong lensing studies) provide us with a unique 
way to probe the statistical properties of the underlying den- 
sity distribution. Traditionally, the study of gravitational 
clustering in the quasi-linear and non-linear regimes was 
done by analyzing galaxy catalogues and were plagued by 
the question of bias associated with the galaxy distribution. 
Now, therefore, one needs a prescription for galaxy bias to 
relate the galaxy distribution to the underlying dark mat- 
ter distribution. Gravitational lensing allows us to bypass 
this problem as it directly probes the mass distribution. Pi- 
oneering works in this direction were done by Blandford et 



al. (1991), Miralda-Escude (1991), and Kaiser (1992), based 
on the earlier work by Gunn (1967). 

A key tool for the understanding of cosmological weak 
lensing has been provided by the development of numeri- 
cal techniques to follow null geodesies in a perturbed model 
universe. Such numerical studies typically employ iV-body 
simulations through which ray-tracing experiments are con- 
ducted (see, e.g., Schneider & Weiss, 1988, Jaroszynski et 
al., 1990, Lee & Paczynski, 1990, Jaroszynski, 1991, Babul 
& Lee, 1991, and Blandford et al., 1991). Building on the ear- 
lier work of Wambsganss et al. (1995 and 1997), detailed nu- 
merical studies were done by Wambsganss, Cen & Ostriker 
(1998). Other recent studies using ray-tracing have been con- 
ducted by Premadi, Martel & Matzner (1998), Van Waer- 
beke, Bernardeau & Mellier (1999) and White & Hu (2000). 
Recently, Couchman, Barber & Thomas (1999) have devel- 
oped a new technique based on computing the full three- 
dimensional shear at locations along the lines of sight. Their 
method combines the computed shear matrices to produce 
the required lensing statistics for sources at the selected red- 
shift. The method is complimentary to ray-tracing and has 
been fully implemented. Excellent agreement between the 
results obtained using Couchman et al.'s (1999) method and 
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analytical predictions have been reported by Barber (2002) 
and Barber & Taylor (2002). 

On the other hand, analytical approaches have also been 
developed to elucidate the link between the statistical distor- 
tion of distant sources by weak lensing and the properties of 
the large-scale structure of the universe. Most of these ana- 
lytical works have focussed on large smoothing angles where 
perturbative calculations are valid (e.g., Villumsen, 1996, 
Bernardeau et al., 1997, Jain & Seljak, 1997, Kaiser, 1998, 
Van Waerbeke, Bernardeau & Mellier, 1999, and Schneider 
et al., 1998). Small angular scales are more difficult to han- 
dle since they probe the non-linear regime where the density 
field can no longer be described through perturbative expan- 
sions which start behaving badly (e.g., Valageas, 2002b). 
Thus, for ongoing surveys with a small sky coverage per- 
turbative results are unlikely to be sufficient. Therefore, we 
need theoretical predictions at small angular scales where 
the signal will dominate the noise. 

Unfortunately, we lack a complete understanding of 
gravitational clustering in the highly non-linear regime since 
it has not been possible so far to solve the Vlasov-Poisson 
system which describes the collisionless gravitational dy- 
namics. Nevertheless, several phenomcnological models have 
been developed to describe the non-linear density field. An 
interesting family of such models, often called "hierarchical 
models", assumes a tree hierarchy for many-body correla- 
tion functions. Then, each model is specified by the way it 
assigns weights to trees of the same order but of different 
topologies (Fry, 1984, Schaeffer, 1984, Bernardeau & Scha- 
effer, 1992, and Szapudi & Szalay, 1993 and 1997). Note that 
the evolution of the two-point correlation function in all such 
approximations is left arbitrary. On the other hand, inde- 
pendent studies by various authors (Hamilton et al., 1991, 
and Peacock & Dodds, 1994 and 1996) suggest an accurate 
fitting formula for the evolution of the two-point correlation 
function based on simple non-local scaling arguments. Com- 
bining a hierarchical Ansatz with such a non-local scaling 
for the two-point correlation one obtains a complete sta- 
tistical description of gravitational clustering in the highly 
non- linear regime (e.g., Valageas & Schaeffer, 1997, Munshi 
et al., 1999a, Munshi, Coles & Melott, 1999b,c, and Munshi, 
Melott & Coles 1999d). 

Such an approach was initiated by Hui (1999) for the 
study of the three point correlation function for convergence 
maps and was later extended to compute the full probabil- 
ity distribution function (Valageas, 2000a, b, and Munshi & 
Jain, 2000, 2001), the associated bias (Munshi, 2000) and 
the cumulant correlators associated with such distributions 
(Munshi & Jain, 2000). See also Munshi & Wang (2003) 
for the extension of these studies to cosmological scenar- 
ios with dark energy. Next, such studies were extended to 
handle the compensated filter appropriate for the aperture- 
mass, M a p, which is more useful in weak lensing surveys 
(Bernardeau & Valageas, 2000). Error estimations were also 
performed based on such a hierarchical Ansatz to optimize 
survey strategies (Munshi & Coles, 2003), while implications 
for Type la Supernova? studies were carried out by Valageas 
(2000a) and Wang, Holz & Munshi (2002). 

So far most analytical calculations have focussed on con- 
vergence maps where analytical results are simplified due to 
the fact that the convergence is a scalar. Similar progress in 
the case of the shear was lacking due to the spinorial na- 



ture of shear fields. However from an observational point of 
view, shear maps are a direct outcome of suveys whereas 
map-making for the convergence field can be problematic 
due to the non-trivial topology of the survey area. There- 
fore, extending previous studies of weak gravitational lensing 
based on hierarchical models, we introduce in this article a 
model, which we also call a "stellar model," to evaluate the 
probability distribution function, PDF, of the shear. We test 
our analytical predictions against numerical simulations and 
we show that our simple model provides a good description 
of the shear PDF at small angular scales. Our approach is 
complimentary to studies related to lower order statistics of 
the shear field (Takada & Jain, 2003, Zaldarriaga & Scocci- 
marro, 2003, Schneider & Lombardi, 2003, and Bernardeau, 
Mellier & van Waerbeke, 2002) which also aim to detect 
non-Gaussianity from ongoing surveys. 

This paper is organized as follows. Section 2 briefly de- 
scribes the weak gravitational lensing effects as measured 
through the convergence, k, and the shear, 7. In Section 
3, we recall the relationship between PDFs and cumulants 
(which are the basic tools of our analytical method) and we 
present our model for the many-body density correlations. 
This allows us to compute in Section 4 the PDFs for the 
shear components and for the shear modulus. Section 5 con- 
tains a brief explanation of the method used to obtain the 
lensing statistics from the numerical simulations. In Section 
6, we make the comparison between these numerical results 
and our anaytical predictions, and we conclude in Section 7 
by discussing our results and future work. 



2 DISTORTIONS INDUCED BY WEAK 
GRAVITATIONAL LENSING 

2.1 Shear tensor 

The gravitational lensing effects produced by density fluctu- 
ations along the trajectory of a photon lead to an apparent 
displacement of the source and to a distortion of the im- 
age. Thus, light coming from a direction 9 is deflected by a 
small angle 88. However, the observable quantities are not 
the displacements 89 themselves but the distortions induced 
by these deflections, which are given by the symmetric shear 
matrix (e.g., Jain, Seljak & White, 2000): 



Here x is the radial comoving coordinate (and Xs corre- 
sponds to the redshift, z s , of the source): 

7T dz 

dx = " (2) 

^Qa + (1 - fi m - 00(1 + z) 2 + n m (l + zf 

while the angular distance T> is defined by: 
■sin* (I l-n m -n A | 1/2 Hox/c) 



H 



i-o. 



(3) 



where sin^ means the hyperbolic sine, sinh, if (1 — Q m — 
Sl\) > 0, or sine if (1 - fi m - Q A ) < 0; if (1 - fi m - Q A ) = 0, 
then T>(x) = X- The gravitational potential, 0, is related 
to the fluctuations of the density contrast, 8, by Poisson's 
equation: 
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H, 



A4> = -n m (1 + z) 5 with <5(x) 



P(x) 



( 4 ) «, = ^ * ^2 J*' dxw( x , Xs )8( X ,V$). 



where p is the mean density of the universe. In eq.(l) we used 
the weak lensing approximation; the derivatives ViVj<p(x) 
of the gravitational potential are computed along the un- 
perturbed trajectory of the photon. This assumes that the 
components of the shear tensor are small but the density 
fluctuations 5 can be large (Kaiser, 1992). The shear tensor 
$ij is usually decomposed into its trace, ft, and the shear 
components 71, 72, defined by: 

$1,1+ $2,2 



and 



71 = 



$ 2 



72 



-$1,2 , 7 = 71 + i 72. 



(5) 



(6) 



Note that since the shear components, 71 and 72, and their 
signs can be exchanged through appropriate rotations of the 
coordinate axes the joint probability distribution function 
(PDF) £"(71,72) obeys: 



P(7i,72) =P(|7i|,M) =P(l72|,l7iD- 



2.2 Convergence: n 



(7) 



Using eq.(l) and eq.(4) one can show (Bernardeau et al., 
1997, and Kaiser, 1998) that the convergence along a given 
line of sight is: 



with: 



Hi V(x)V(xs-x) 
c 2 V( Xs ) 



(1 + *) 



(8) 



(9) 



where z corresponds to the radial distance X- Thus the con- 
vergence, k, can be expressed very simply as a function of 
the density field; it is merely an average of the local density 
contrast along the line of sight. Then, we can see from eq.(8), 
that there is a minimum value, Kmin(z s ), for the convergence 
of a source located at redshift z s , which corresponds to an 
"empty" beam between the source and the observer (S = — 1 
everywhere along the line of sight): 



3n n 



dx w(x,Xs)- 



(10) 



Following Valageas (2000a, b) it is convenient to define the 
"normalized" convergence, k, by: 



dx w 8, with w — 



w(x,Xs 



Jo 3 d X w(x,Xs) 



(11) 



which obeys k > — 1. Here we introduced the "normalized 
selection function," w( X ,Xs)- One interest of working with 
normalized quantities like k is that most of the cosmological 
dependence (on fi m ,f2A and z s ) and the projection effects 
are encapsulated within K m in, while the statistics of k (e.g., 
its PDF) mainly probe the deviations from Gaussianity of 
the density field which arise from the non-linear dynamics 
of gravitational clustering. If one smoothes the observations 
with a top-hat window in real space of small angular radius, 
9 S , one rather considers the filtered normalized convergence 
k s (where the subscript "s" refers to "smoothed"): 



(12) 



Here i9 is a vector in the plane perpendicular to the line of 
sight (we restrict ourselves to small angular windows) over 
which we integrate within the disk t?j < 6 S ; we note this by 
the short notation J^ 3 . Thus x is the radial coordinate, while 
T>-d is the two-dimensional vector of transverse coordinates. 
Eq.(12) clearly shows that the convergence k s is actually 
an average of the density contrast over the cone of angular 
radius 6 S . 

In the following, it will be convenient to work in Fourier 
space. Thus, we define the Fourier transform of the density 
contrast by: 



*(x) 



/ 



dk e lk x <5(k) 



(13) 



where x and k are comoving coordinates. Then, eq.(12) also 
reads: 

*•= j d X Hx,Xs) j dke ik « x W(k ± V6 s ) J(k), (14) 

where ku is the component of k parallel to the line of sight 
and kx is the two-dimensional vector formed by the compo- 
nents of k perpendicular to the line of sight. Here we intro- 
duced the Fourier form W(k±V6 s ) of the real-space top-hat 
filter of angular radius 8 S : 



w(k ± ve s ) = 



7V9 2 S 



k ± V6 s 



(15) 



where Ji is the Bessel function of the first kind of order 1. 
If we choose another filter (e.g., a Gaussian window rather 
than a top-hat) the expression (14) remains valid and we 
simply need to use the relevant Fourier window W(k±V9 s ). 

2.3 Shear: 7 

For the shear, 7, defined in eq.(6), we can obtain expressions 
which are similar to eqs.(12)-(14). Thus, from eq.(l) and 
cq.(6) we obtain in Fourier space for the normalized shear 

7 = 7/l«min|: 



7 : 



dxw(x,Xs) / dke 



ik.x k\ — k\ + 2ikiki 



kl + k 2 2 



*(k) (16) 



where ki and ki are the components of k along the two 
orthogonal axes perpendicular to the line of sight which also 
define the directions (1, 2) used in Subsection 2.1 for the 
components Qij of the shear tensor. In eq.(16) we used again 
the small-angle approximation. Of course, for the smoothed 
shear y s we get, in a similar fashion to eq.(14): 



/ dx» / dke' 
xW(k ± V6 s ) <5(k). 



k\ - k\ + 2ikik 2 



(17) 

Going back to real-space, we also obtain from eq.(17) for a 
top-hat filter the expression: 

%=- r dx w r f r ^ e - ( i 8) 

v " & s " 

where 1? is a vector in the plane perpendicular to the line of 
sight of length t? and which makes the polar angle, a, with 



© 0000 RAS, MNRAS 000, 000-000 



4 P. Valageas et al. 



the 1-axis (the 2-axis has a = 7r/2). Note that we integrate 
over i? over all space outside of the disc of radius 6 S . Thus, 
the comparison with eq.(12) shows that while the smoothed 
convergence only depends on the matter within the cone 
formed by the angular window, 6 S , the smoothed shear only 
depends on the matter outside this cone. 



3 CUMULANTS, PDFS AND MODELS FOR 
THE MANY-BODY DENSITY 
CORRELATIONS 

In this paper, we wish to evaluate the probability distribu- 
tion function, PDF, of the shear, 7. To this order, following 
the approach developed in Valageas (2000a, b) and Munshi 
& Jain (2000, 2001), we shall first compute the moments (or 
more precisely the cumulants) of the shear and then obtain 
the associated PDF. Therefore, we first recall the standard 
relationship between the PDF and the moment and cumu- 
lant generating functions. 



3.1 Generating functions 

The statistical properties of a random variable s can be ob- 
tained from its moment or cumulant generating functions, 
which are widely used in statistics (see also Balian & Scha- 
effer, 1989). Thus, the PDF V{s) can be derived from the 
generating function, ip(y), defined from the moments, (s p ), 
(provided they are finite) by: 



E ^ ^ y p 

p=0 



p< 



since we have the inverse Laplace transform: 



V{S) 



dy_ 

2m 



(19) 



(20) 



Indeed, the generating function, tp(y), is simply the Laplace 
transform of the PDF V{s): 



= 1 



As e- sv V(s). 



(21) 



One can check that expanding the exponential in eq.(21) re- 
produces eq.(19). If the moments, (s p ), diverge one can still 
define tp(y) from eq.(21) but the expansion at y = of the 
generating function becomes singular. In practice, one usu- 
ally introduces the generating function, 3>(y), defined from 
the cumulants, {s p ) c , by: 



My) 



.i)p 



P =i 



pi 



(s p ) c y p = \n[^{y)]. 



(22) 



As is well-know, $(y) is simply the logarithm of tj)(y). For 
our purposes, we shall actually use the normalized generat- 
ing function, ip(y), defined by: 



<p(v) = E 



(-If" 1 {S^c 



pi 



p-i 



y p with £ s = (s 2 ) c . 



(23) 



The comparison of eq.(23) with eq.(22) yields ip(y) = 
exp[— <p(yt;s)/£,s], so that eq.(20) now reads: 



V{s) 



-io 



Ay 



\sy-<p(y)]/is 



(24) 



As seen in Section 2, the convergence, k, or the shear, 7, 
are given by a linear integral along the line of sight over the 
density field. In other words, they can be seen as the super- 
position of independent layers which each contribute to k or 
7. Indeed, within the small-angle approximation the corre- 
lations of the density field will only appear within the trans- 
verse directions. Then, the direct calculation of the PDF 
of such a sum over the redshift of the lenses would yield 
an infinite number of convolution products which makes 
it intractable. By contrast, the cumulant generating func- 
tions simply add when different layers are superposed (since 
Laplace transforms change convolutions into ordinary prod- 
ucts). This property makes the generating function ip(y) a 
convenient tool to deal with projection effects along the line 
of sight. This remark is the basis of the method introduced 
in Valageas (2000a, b). Moreover, the expansion (23) pro- 
vides a simple way to compute <p(y) from the moments (or 
cumulants) of the random variable one is interested in. This 
approach can be used for any quantity which is linear over 
the density field, like the convergence, k, studied in Valageas 
(2000a, b), the aperture-mass, M ap , investigated in details 
in Bernadeau & Valageas (2000) or the shear, 7, we focus 
on in this paper. 

3.2 PDF for the density field 

It is clear that in order to derive the PDF of the shear, 7, 
we need the properties of the underlying density field. A 
first measure of the statistical properties of the density field 
is given by the PDF V(5r) of the density contrast within 
spherical cells of radius R and volume V: 



dx 



(25) 



Note that the cumulants, {S R ) C , can also be written in 
terms of the many-body connected correlation functions, 
£ p (xi, .., x p ), as (e.g., Peebles, 1980): 



P>2: (^>c = Ip = j^^^6.(xi,..,x J> ) 1 
with: 

£ P (xi, .., x p ) = (<$(xi)..<S(x p )) c . 



(26) 



(27) 



Next, we define the associated normalized generating func- 
tion, ip(y), introduced in eq.(23) and its coefficients, S p , by: 



<p(v) = E 



(-1) 



p-i 



P 



S P y p 



(28) 



(29) 



which also yields V(Sr) through eq.(24): 

/ + ioo j _ 
d V_ e [SRV-<P(v)]/(2 m 
■ioo 27r *£ 2 

Since the convergence, k a , is the average of the density con- 
trast, 5, with a smooth positive weight over the relevant 
cone, its cumulants can be estimated from those of the den- 
sity field in a robust and simple way, by neglecting their ge- 
ometrical dependence. Therefore, one merely uses the prop- 
erties of the density field, 5r, smoothed over spherical cells 
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and does not take into account the detailed angular depen- 
dence of the many-body correlation functions. As shown in 
Valageas (2000b), this procedure is indeed sufficient to re- 
cover the properties of the convergence. In fact, up to a good 
accuracy, one may even express the PDF, V(k 3 ), directly in 
terms of the PDF, V(Sr), of the density contrast, by us- 
ing the approximation <fia s {y) — <p{y) (see Valageas 2000b). 
Here tp (also defined in eq.(28)) is the normalized generating 
function associated with the density contrast Sr at the scale 
and redshift which contribute most to the signal, whereas 
ifft s is the normalized generating function associated with 
the normalized convergence, k 3 , as defined by eq.(23). This 
approximation clearly shows that the "normalized conver- 
gence," k s , mainly probes the non-Gaussianities of the den- 
sity field (i.e., its very structure) brought by the non-linear 
gravitational dynamics, while the overall amplitude of the 
fluctuations (which gives the dependence on the cosmolog- 
ical parameters and the source redshift) is encoded within 
Kmin, as we announced earlier, below eq.(ll). 



3.3 Models for the many-body density 
correlations 

For more intricate filters, this simple approximation is no 
longer sufficient. This is clearly the case for the shear, 7, 
since eq.(18) shows that the relevant volume is now the ex- 
terior of the cone of radius 9 S , and has a geometry which is 
significantly different from spherical cells. More importantly, 
the factor e l2a yields a weight of mean zero and which is both 
positive and negative and varies with the direction in the 
transverse plane. This implies that one cannot evaluate the 
cumulants of 7 from these moments of the density field over 
spherical cells; the natural estimate would be zero, since the 
angular mean of the weight, e t2a , vanishes. The same prob- 
lem appears for the study of the aperture mass, M ap , which 
also involves a compensated filter. Therefore, we need an ex- 
plicit model for the angular dependence of the many-body 
correlations. 

As shown in Bernardeau & Valageas (2000), in the case 
of the aperture mass, M ap , good results can be obtained at 
small angles by using a "minimal tree-model" for the con- 
nected correlations of the density field. Let us first recall that 
a general "tree-model" is defined by the property (Schaeffer, 
1984, and Groth & Peebles, 1977): 



& (xr , . . , x p ) = Q £ a) E II & < Xi ' x * ) 



(30) 



(«) 



where (a) is a particular tree-topology connecting the p 
points without making any loop, Q p a ^ is a parameter as- 
sociated with the order of the correlations and the topology 
involved, t a is a particular labeling of the topology, (a), 
and the product is made over the (p — 1) links between the 
p points with two-body correlation functions. We show in 
Fig.l the three topologies which appear within this frame- 
work for the 5-point connected correlation function. 



3.3.1 Minimal tree model 

A peculiar case of the tree-models described by (30) is the 
minimal tree-model where the weights Q p °^ are given by 
(Bernardeau & Schaeffer, 1992): 
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(a) (b) (c) 





Figure 1. The various topologies one can build for the 5-point 
connected correlation function within the framework of a tree- 
model as in (30). Graph (a) is a "stellar diagram" (i.e., (p — 1) 
points arc linked to a central point), while graph (b) is a "snake 
diagram" (i.e., one successively visits all points along one curve 
with only two end-points) and graph (c) is a "hybrid diagram" . 



Qp 



(a) 



n 



(31) 



vertices of (a) 



where v q is a constant weight associated to a vertex of the 
tree topology with q outgoing lines. Then, for an arbitrary 
real-space filter, -F(x), which defines the random variable s 
as: 



dxF(x)<5(x) and £ s = (s 2 ), 



(32) 



it is possible to obtain a simple implicit expression for the 
normalized generating function, ip s (y) (see Bernardeau & 
Schaeffer, 1992, and Jannink & Des Cloiseaux, 1987): 



r(x) = 



V / dx F(x) 



I 



C[r(x) 



t(x)C'[t(x)] 



dx'F(x')6(x,x') C'[r(x')] 



(33) 



(34) 



where the function £(t) is defined as the generating function 



cw=E 



(-I)* 



with 



1. 



(35) 



Of course, the generating function, <p s (y), depends on the 
filter, F, which defines the variable, s. For the smoothed 
density contrast, Sr, defined in eq.(25), the real-space fil- 
ter, F(x), is the top-hat of radius R normalized to unity. In 
this case, a simple "mean field" approximation which pro- 
vides very good results (Bernardeau & Schaeffer, 1992) is 
to integrate r(x) over the volume V in eq.(34) and then to 
approximate r(x) by a constant r. This leads to the simple 
system: 



<p(y) 



y 



C(r) 



r C'(r) 



(36) 
(37) 



Here and in the following we simply note the generating 
function associated with the density contrast Sr as (p(y) 
and we retain suffixes for the generating functions associated 
with other quantities (e.g., Lpx lls for the shear component, 
7is)- 

In the case of the aperture mass, M ap , one can still 
perform the resummation (33)-(34), if the underlying den- 
sity field is described by such a minimal tree-model, but one 
can no longer use the "mean field" approximation (36)-(37) 
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because the 2-d compensated filter F(-&) now strongly de- 
pends on the polar distance |$| (in the 2-d transverse plane). 
Such a study is described in detail in Bernardeau & Valageas 
(2000) where it is shown that this minimal tree-model pro- 
vides good results when compared with numerical simula- 
tions. 

Of course, we can apply the same procedure, developed 
in Bernardeau & Valageas (2000) for M ap , to the shear, 
7. We simply need to use the corresponding real-space fil- 
ter, F. From eq.(18) we see that this is F(S) = -H(d > 
6 S ) cos(2a)/7n9 2 for 71 and F(d) = -H(0 > 9 S ) sin(2a)/7n9 2 
for 72, where H is the Heaviside function with obvious 
notations. Then, the resummation (33)-(34) yields a two- 
dimensional non-linear implicit equation which is not very 
convenient for numerical purposes. Therefore, it is interest- 
ing in this case to investigate whether a simpler model for 
the many-body correlations of the density field could still 
provide good predictions for the PDF, V(-yi), of the shear 
component, 71. 

3.3.2 Stellar model 

Thus, we introduce in this paper a "stellar-model" , which is 
a particular case of the tree-models defined in (30), where we 
only keep the stellar diagrams (e.g., the graph (a) in Fig.l for 
the 5-point connected correlation). Thus, the p— point con- 
nected correlation £ p of the density field can now be written 
as: 



fp(xi,..,X p ) = ^ Jj£ 2 (Xi,Xj). 

i=l jjti 



(38) 



The advantage of the stellar-model (38) is that it leads to 
very simple calculations in Fourier space. Thus, let us define 
the power-spectrum, P(k), of the density contrast by: 



(5(ki)<5(k 2 ))=fo(ki+k 2 ) P(ki), 
where So is Dirac's distribution, whence: 



dk e lk x P(k). 



(39) 



(40) 



Then, eq.(38) reads in Fourier space: 

~ p 
(5(k 1 )..8(k p )} c = ^ fo(kr + .. + k p ) /'!/.-. :. (41) 

i — l jj^i 

Of course, the Dirac factor, <5o(ki + .. + k p ), simply trans- 
lates the fact that the many-body correlations are invariant 
through translations. Therefore, in the following we shall 
use eq.(41) for the connected correlations of the density 
field which appear in the calculation of the cumulants of 
the shear, 7. 



3.4 Amplitude of the density correlations 

The stellar-model introduced in Section 3.3.2 provides the 
angular dependence of the many-body correlations of the 
density field. However, we still need to determine their over- 
all amplitude, that is the non-linear power-spectrum, P(k), 
and the coefficients, S p , which enter eqs.(38)-(41). First, we 
obtain the non-linear evolution of the power-spectrum from 
the prescription given by Peacock & Dodds (1996), which is 



also a good fit to numerical simulations. Second, we deter- 
mine the coefficients, S p , as follows. We characterize the am- 
plitude of the density correlations at scale 7? and redshift z 
through the parameters S p defined in eq.(28), which give the 
amplitude of the many-body correlations over spherical cells 
of radius R. As seen from eq.(28)-(29) they also define the 
PDF, V(Sr). We could obtain these coefficients, S p , through 
a parameterization of their generating function, <p(y). How- 
ever, we prefer to specify (p(y) through a function f (t) which 
is defined by the implicit system (36)- (37). Note that this 
procedure is merely a convenient parameterization of <p(y) 
and it does not assume by itself a minimal tree-model. The 
results recalled in Section 3.3.1 simply mean that in case the 
actual many-body correlations exactly obey a minimal tree- 
model, the function £(t) defined from eqs.(36)-(37) could 
also be interpreted as a good approximation to the (slightly 
different) function C„(t) defined in eq.(35) as the generating 
function of the vertices v p . 

The reason we prefer to parameterize the coefficients, 
S p , through £(t) rather than directly through tp(y) is that 
this procedure was already checked to provide good results 
in previous works. Thus, it was seen in Bernardeau & Schaef- 
fer (1992) to give good predictions for the galaxy and matter 
correlations, while in Bernardeau & Valageas (2000) it was 
shown to yield good predictions for the PDF, 7>(M ap ), of 
the aperture- mass. Since the aperture-mass is closely related 
to the shear, 7, we can expect this procedure to give good 
results for the PDF, "P(7), too. Note however that the calcu- 
lations developed in Bernardeau & Valageas (2000) used the 
minimal tree-model recalled in Section 3.3.1 rather than the 
stellar-model of Section 3.3.2 for the angular dependence of 
the many-body correlations. Nevertheless, we checked nu- 
merically that both models actually give close results for 
■P(M ap ). Thus, as in Bernardeau & Valageas (2000) we pa- 
rameterize Q{t) as: 



CM 



1 + 



i, 



(42) 



where we have kept the usual notation, k (not to be confused 
with the convergence), for the free parameter which enters 
the definition of C( r )- 

Next, we choose the value of this parameter, k, so as 
to reproduce the skewness, S3, in both the highly non- linear 
and quasi-linear regimes, using the relation: 

* = ^h- <«> 

For S3 > 3 we have k > while for < S3 < 3 we have 
k < —1. In the limit S3 — > 3 the function £(t) defined in 
eq.(42) goes to f(r) — > e~ T — 1 (and \k\ — » 00). We take 
for the skewness in the highly non-linear regime, S^ 1 ", the 
prediction of HEPT (Scoccimarro & Frieman, 1999) and for 
the quasi-linear regime, S® h , the exact result obtained from 
perturbative theory: 

4-2" „ot, 34 
7 



oQL _ 

°3 — 



(n + 3). 



(44) 



1 + 2 n+1 ' 

Here n is the local slope of the linear power-spectrum at 
the typical wavenumber k s probed by the observations at 
redshift z. From eq.(14) and eq.(17), we define this typical 
wavenumber as: 
1 



k s (z) = 



V(z)9 s 



(45) 
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Finally, we introduce the power, A 2 (k), per logarithmic in- 
terval defined as: 

A 2 (k,z) =4ivk 3 P(k,z). (46) 

Then, for intermediate regimes defined by 1 < A 2 (k s , z) < 
Avir(z), we use the simple linear interpolation: 

Ss(z) = + A ^Si ( 5 ? L - S^) . (47) 

For A 2 (k s , z) < 1 (i.e. the quasi-linear regime), we take S3 = 
S^ L , while for A 2 (fc s ,z) > A v i r (z) (i.e., the highly non- 
linear regime), we take S3 = S^ L . Here A v i r is the density 
contrast at virialization given by the usual spherical collapse 
(thus A vir ~ 178 for a critical- density universe). Indeed, the 
threshold A 2 > A v i r describes the highly non-linear regime 
where most of the matter at scale 1/k has collapsed into 
non-linear structures. 

Finally, as can be seen from eq.(41), the coefficients S p 
and S p are related by: 



5p — Sp J 
Jo 



dt 3r 



Jf A 2 (k)F(kRr-^l 
Jf A 2 (k)F 2 {kR) 



p-i 



(48) 



where we introduced the Fourier transform, F(kR), of a 3-d 
top-hat of radius R: 

F(kR) =jf e lk x = 3 Sin{kR) ~^ C ° S(kR) - (49) 

However, in the following we shall use the simple approxi- 
mation: 

S P ^ Sp. (50) 

Alternatively, we may define the function <p(y) obtained 
from (36)- (37) and our choice of £(t) as the generating func- 
tion of the coefficients S p , rather than S p , through its Taylor 
expansion at y = 0. 

It is interesting to note that the implicit system (36)- 
(37) appears naturally in the quasi-linear regime. Indeed, as 
shown in Bernardeau (1994) and Valageas (2002a), in this 
regime the generating function <p{y) exactly obeys the rela- 
tion (36)-(37), where £(t) is closely related to the spherical 
dynamics for the non-linear density contrast. However, the 
latter is of the form (42) only in the limit Q m — > with no 
smoothing (then, the parameter k is equal to 3/2). Never- 
theless, the form (42) still provides a good approximation 
for more general cases where n is adjusted so as to recover 
the exact skewness. Here we note that the system (36)-(37) 
usually yields a branch cut along the negative real axis for 
y < y s with y a ~ —0.1 for ip(y); see Bernardeau & Schaeffer 
(1992). For the function £(t) given in eq.(42) one gets (for 
k > or k < -2, that is S3 > 3/2): 

y ' = -7+2(1^) and T ° = -^T2- (51) 

Such a singularity leads to an exponential cutoff for the 
PDF, V(5 R ) ~ e VsSR/ ^ 2 at large densities. As shown in 
Valageas (2002a), in the case of the quasi-linear regime this 
is actually an artefact and one must go through y s up to 
a second branch for ip(y), which yields a large density cut- 
off for V(Sr) which is shallower than a simple exponential. 
However, it was seen by numerical computations in Valageas 
(2002a) that this feature does not change the shape of the 
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Figure 2. Probability distribution function, 7 3 (7i), for the 
LCDM cosmology, plotted as a function of 71 . The smoothing an- 
gle, 8 a , varies from 1' (upper left panel) to 8' (lower right panel). 
The data-points are the results from the numerical simulations 
and the solid lines correspond to the analytical predictions (62) 
and (64) based on our stellar-model. The dashed lines show the 
Gaussian PDF with the same variance. The source redshift is 
fixed at unity in each case. 



PDF by much in the range of interest. Moreover, in this pa- 
per we are mostly interested in the non-linear regime (which 
largely dominates for 9 ^ 4' and z s = 1) where the singular- 
ity y s is meaningful and the PDF shows indeed an exponen- 
tial cutoff (since it is actually defined from the coefficients 
S p ). 

Finally, we must point out that the simple model de- 
veloped above has no new free parameter. Indeed, the power- 
spectrum evolution is obtained from Peacock & Dodds (1996), 
while all many-body correlations are defined by the stellar- 
model (38) with their amplitude given by the sole coeffi- 
cient 53 (i.e. the skewness of V(5r)) which is exact in the 
quasi- linear regime and obtained from HEPT (Scoccimarro 
& Frieman, 1999) in the non-linear regime. In the form given 
here, our model is obviously consistent with the "stable- 
clustering Ansatz" (e.g., Peebles 1980). However, it could 
incorporate possible deviations from the stable-clustering 
Ansatz through the non- linear power-spectrum, P(k), or 
through an additional dependence on redshift and scale for 
the skewness, ^3, in the non- linear regime. 



4 COMPUTATION OF THE PDFS 

4.1 PDF of the shear components, 71 and 72 

We can now use the model developed in Section 3 to describe 
the density field in order to compute the PDF, P(7i), of the 
shear component, 71. First, we obtain from eq.(17): 
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Figure 3. As for the previous Figure, but for the OCDM cos- 
mology. 



71.- 



= J dx w J < 



dk e ifc " x cos 20 W(k±Ve s ) 5(k), (52) 



where <f> is the polar angle of the transverse vector kx- For 
the component 72 we would have sin 2cj> in place of cos 2<f>. 
Then, the cumulant of order p reads: 



<7f s )c 



rxs p r p 

= / Il dXi {bi I II dkj cos2( t>i W(k ±j V9, 
Jo i=1 J j=1 



X . f[e 

\i=i / 

Since the correlation length (beyond which the density cor- 
relations are negligible) is much smaller than the Hubble 
scale, c/H(z) (where H(z) is the Hubble constant at red- 
shift 2), we can simplify eq.(53) as: 

pXs roo P p P 

<7f s )c= / dxitSf / l[dxi / YldkjcoBfyjWfajVe, 

Jft J —00 j_2 j = 1 

x ^Qe'*!!'*^ 5 p fc(k 1 + ..+k p )P(fc 2 )..P(fc p ),(54) 



, ik «lXl 



(<y(ki)..<y(k p )> c 



) 

(53) 



where we used eq.(41). Here the coefficients S'p(^i) are eval- 
uated at the typical wavenumber, k s (zi), defined in eq.(45), 
probed by the smoothed shear at the redshift z\ along the 
line of sight. Using the usual small-angle approximation 
(i.e. P(kj) ~ P(k±j)), we can perform the integration over 
X2, -,Xp and --.fcllp, which yields: 

p 



<7L)c 



2tt 



(27TW)) 



3=1 



x S* p fo(k±i + .. + k ±P )P(fc±2)..P(fcxp). (55) 
Next, we write the Dirac factor as: 



fo(kxi + .. + kx p ) 



dt 

(2tt)' 



,»t.(k_ L1 + ..+k Lp ) 



(56) 



so that eq.(55) reads: 
{Y ls )c = f^S p (2nw) p 



dt 



dki..dk p 



x J| cos 2^ WikjVOsy^ 00 *^-^ x P(fe 2 )..P(fe p )(57) 
2=1 

where we have dropped the subscript, _L, for the 2-d trans- 
verse wavenumbers kx and we have noted (f> the polar angle 
of the vector t. The integration over the angles <f>j yields: 



itkj cos(<fij ~4>) 



d(f>j cos 2(f>j e 



so that eq.(57) becomes: 



= -27Tcos20 J 2 {tkj) (58) 



<7f s )c 



dx S p w 



■ r 1 dtt(ve s ) 

Jo 



" /"g(-cos20r 



f 

Jo 



dfci k 1 w(k 1 ve s )j 2 (tk 1 ve s ) L fl (t,x) r 



(59) 



where we have made the change of variable t — ► tV6 and we 
have introduced J 71 (t, x) defined by: 



J 71 (*, X) 



4ty' 



poo 

/ dk kw(kve a )j 2 (tkve s )p(k) 

Jo 

A 2 (k,z) 



[°° dk 

71 L T 



w{kve 3 )j 2 {tkve a ). (60) 



For the angular top-hat filter, W , given by eq.(15) in Fourier 
space, we can perform the integration over fci in eq.(59), 
using the standard properties of Bessel functions, and we 
obtain: 



<7f s >c = 



dx r t / " — S> (»cos20) p / 7 7 1 . (61) 
Ji 1 Jo n 

Here we have dropped the minus sign in front of the factor 
cos 2(j> because all odd moments vanish, as seen by integra- 
tion over <j>. Thus, we recover the fact that the PDF, "P(7i), 
is even, which was obvious from the definition of the shear 
since there is no preferred direction in the system. 

Finally, the expression (61) allows us to obtain the gen- 
crating function <fij ls (y) of the shear component 71 s , defined 
as in eq.(23), using the resummation (28) (where we replace 
S p by S p ; see eq.(50)): 

/■2tt 

; ; /it 1 

<piis (v) 



(62) 



Jo Ji 1 Jo 71 



X ifi [ y cos 2(f> w; z 

S7ls 



Here we have introduced the variance £ 71s = (7L). From 
eq.(61) it is given by: 

£ 71s = / d X / —J -f (™cos20) 2 / 71 

1 f Xs , .2 f 00 dk A 2 (k,z) , 2 
= 2 J dxW 77 J 1 k — W ( kV6 *) ■ ( 63 ) 

Here we have used S 2 = 1; see eq.(41), which is still verified 
by our parameterization (42). Of course, eq.(63) is exact, 
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within the small-angle approximation, and the last equality 
in eq.(63) shows that we recover the property (7i s ) = («?} /2- 
For the shear component, 72, we obtain the same equations 
with cos 2(f) replaced by sin 2(f). This leaves the cumulants 
(ifs )c = (7is)c unchanged and we get the same generating 
function. Therefore, we recover P(~/i) = ^(72)- Finally, the 
PDF is given by the inverse Laplace transform (24): 



/-He 
-ioo 



dy 



2ivi^ 1 



= [71 



,(y)]/^ 



(64) 



The generating function <p(y) obtained in Section 3.4 shows 
a branch cut along the negative real axis for y < y s ; see 
eq.(51). Then, we see from eq.(62) that the generating func- 
tion (fy ls (y) for the shear component 71 shows two sym- 
metric branch cuts along the real axis, for y < —y s ,j l3 and 



with: 



Vs 



— mm 

z.t 



7 71 m) 



(65) 



Let us recall here that the generating function tp(y; z) explic- 
itly depends on the redshift z along the line of sight through 
the skewness Si,{,z) of the density contrast at the typical 
wavenumber k a (z) probed at this redshift, see Section 3.4. 
Hence its singularity y s (z) also depends on the redshift z. 
Next, the singularities of <^> 71s (y) at ±J/ S , 71s also imply an 
exponential tail for V^u) at large |7i s |: 



| 7 i s |>1: P(7i s )~e 



-!/s,7 ls l71sl/S 7 



(66) 



which translates into similar exponential tails for the shear 

710 = |Kmin|7ls- 

4.2 PDF of the shear modulus, |7| 

From the smoothed shear components, 7i 3 and 72s, we may 
define the smoothed modulus, 7 s j, as: 



\%\ = \lis + iT2s\ = \J 7? 3 + 72 



(67) 



Note that eq.(67) means that we take the modulus of the 
shear after smoothing of the two components, 71 and 72. 
Since smoothing and taking the modulus do not commute, 
the quantity |7,| defined in eq.(67) is not the smoothed shear 
obtained by first taking the modulus I7I = I71 +172! at each 
point on the sky and second applying a top-hat smoothing to 
this map. In this paper we shall only consider the smoothed 
modulus |7 S | as defined in eq.(67). Indeed, since taking the 
modulus is a non-linear process one cannot directly apply 
the method used in this paper to derive the PDF of the 
modulus 7 defined by first taking the modulus of the shear. 
By contrast, we can easily study the smoothed modulus |7 S | 
defined in eq.(67) because it can be directly obtained from 
the properties of the smoothed shear components, 7i s and 
72s- Indeed, from eq.(67) we obtain the relation: 



V{\ls\) 



r 

Jo 



da V(\i a \ cosa, \y s \ sin a), 



(68) 



where we have introduced the joint PDF, V{yi s , 72s), of the 
two shear components, 7i s and 72s. Of course, following the 
method used in Section 4.1 for the shear components we 
can obtain the joint PDF, ^(tls, 72s), from the cumulants, 
(7is7 2 ' s ! )c- Thus, proceeding along the lines of Section 4.1 
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Figure 4. Probability distribution function, 7-'(|7|), for the 
LCDM cosmology, plotted as a function of |7|. The smoothing 
angle, 8 S , varies from 1' (upper left panel) to 8' (lower right 
panel). The data-points are the results from the numerical simu- 
lations and the solid lines correspond to the analytical predictions 
(68), (74) and (75) based on our stellar-model. The dotted-dashed 
curve assumes that the two shear components, 71 and 72 , are un- 
corrected but their PDF is given by our non-linear prediction; 
this corresponds to cqs.(76) and (77). The dashed lines are the 
predictions obtained for a Gaussian density field with the correct 
non-linear variance (this yields a Raylcigh distribution for | — y | ) . 
The source redshift is fixed at unity in each case. 



the analog of eq.(61) now reads: 

<tSt£>c = p d X § P1+P2 J°° I p ^ 

x (-cos2<?(>) pl (-sin20) P2 Z^ 2 " 1 . (69) 

Then, the generating function, &j ls ,*, 2s (yi, 2/2), defined as in 
eq.(22) by: 



= E ( X n 7 (C 1 ^^ 2 (70) 

z — ' Pl!p2! 

P1,P2=0 

can be written from eq.(69) as: 

xV-S P M P (yi cos 20 + y 2 sin 2<f>) p . (71) 
^ — ' p\ 
P =i 

Using the resummation (28) we get: 

x ip (— (yi cos 2<f> + 2/2 sin20)7 71 w; z) . (72) 

Next, from the cumulant generating function, $7 ls ,7 2s , we 
define the normalized generating function, ^713,7231 by: 



P1+V2 
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¥>71s,72s (3/1.3/2) = -Cruris ,72s 

and eq.(72) yields: 



yi 3/2 

(.lis (.lis 



(73) 



¥>71 3 ,72„ (3/1.3/2) 



x (p ( (j/i cos 2</> + 2/2 sin 2(£) — 2- to; 2 

S71s 



Then, the joint PDF, P (71s, 72s), is given by: 
dyidj/2 



(74) 



V (713,72s) 



x e 



(2^ 71 J 2 

sJ/1+72s!/2-Vt. 1s 



,(yi.»2)]/4 71 



(75) 



It is interesting to compare our results (74) and (75) with the 
prediction we would obtain if we had assumed that 71 s and 
72s were uncorrelated, but still described by the PDF de- 
rived in Section 4.1. Then, the normalized generating func- 
tion, <Py le ,y 2e (yi,V2) would simply be given by: 



uncorr / \ 
¥ , 71 3 ,72 3 (3/1. 3/2) 



^71a(3/l) + ¥>71 S (J/2). 



(76) 



Substituting this into eq.(75), we obviously recover the prop- 
erty: 



P uncorr (7is,72 S )=P(7is)P(72s). 



(77) 



Therefore, the fact that tp*y le ,y 2 s obtained in eq.(74) cannot 
be written as the sum (76) shows that within our stellar 
model of Section 3.3.2 for the density field, both shear com- 
ponents, 7i s and 72s, are correlated, which is actually quite 
natural. Moreover, the comparison of (75) with (77) will 
allow us to study numerically the influence of these correla- 
tions onto the PDF of the smoothed modulus \%\. 

Finally, we can note that for a Gaussian density field 
we simply have ip(y) = —y 2 /2, see eqs.(26)-(28). This also 
corresponds to S2 = 1 and S p = for p > 3. Hence all our 
results are exact for the case of a Gaussian density field, 
where we simply need to use for the generating function: 
ip Gllues {y) = -y 2 /2. Then, we recover from eq.(62) that the 
shear components 71 and 72 are also Gaussian random vari- 
ables, since we find (^?* uss (j/) oc y 2 . This was actually ob- 
vious from the definition of the shear as a linear integral 
over the density field, see eq.(16). On the other hand, wc 
can see from eq.(74) that in the Gaussian case we have 
^2 S (2/i.3/2) = ^""(yi) + ^""(io), since the cross- 
correlation term yi cos 2<j) x j/2 sin 2(j> vanishes after integra- 
tion over <f>. Thus, we recover the fact that for a Gaussian 
density field the two shear components 71 and 72 along the 
directions (1,2) are uncorrelated. This is not surprising since 
for a Gaussian density field the Fourier modes <5(k) arc un- 
correlated. This yields a Rayleigh distribution for the shear 
modulus |7 S I . 



5 LENSING STATISTICS FROM NUMERICAL 
SIMULATIONS 

We compare in the following Section our analytical results 
for the shear with weak lensing statistics obtained in cosmo- 
logical iV-body simulations. Couchman, Barber & Thomas 
(1999) developed an algorithm for computing the shear in 
three dimensions at locations within the simulation volumes 
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Figure 5. As for the previous Figure, but for the OCDM cos- 
mology. 



Table 1. Cosmological and simulation parameters characterizing 
the different models 





r 




n A 






^survey 


LCDM 
OCDM 


0.25 
0.25 


0.3 
0.3 


0.7 
0.0 


1.22 
1.06 


0'.34 
0'.37 


2°. 6 
2°. 8 



and we have applied this algorithm to the simulations of the 
Hydra Consortium* produced using the 'Hydra' iV-body hy- 
drodynamics code (Couchman, Thomas & Pearce, 1995). 

We used two simulations, a flat cosmology with a cos- 
mological constant and an open cosmology with zero cosmo- 
logical constant. These will be referred to as the LCDM and 
OCDM cosmologies, respectively. Both contained dark mat- 
ter particles only of mass 1.29 x lO 11 /*," 1 solar masses, where 
h expresses the value of the Hubble parameter in units of 
100 km s _1 Mpc -1 . The number of particles in each cos- 
mology was 86 3 . We used a variable particle softening in the 
code with a minimum value in box units of 0.0007(1 + z), 
where z is the redshift of the simulation volume. The simula- 
tion volumes had comoving side-dimensions of 100/i _1 Mpc. 
To avoid obvious structure correlations between adjacent 
boxes, each was arbitrarily translated, rotated (by multi- 
ples of 90°) and reflected about each coordinate axis, and in 
addition, each complete run was performed 10 times in each 
cosmology. The cosmological and simulation parameters of 
both the LCDM and OCDM simulations are given in Table 
1. 

The general procedure for specifying the coordinates 



(http: / /hydra.mcmaster.ca/hydra/indcx.html) 
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for the lines-of-sight and the locations within the simula- 
tions for the computations of the 3-d shear is as described 
by Barber (2002). In the present work, a total of 455 x 455 
lines of sight were used which completely filled the simula- 
tion box immediately preceding the one at redshift 1 and 
allowed regular sampling of the field of view. In addition, 
300 regularly-spaced evaluation locations for the shear along 
each line of sight were specified in each simulation volume 
to adequately sample the varying gravitational potential. In 
the LCDM cosmology, the full field of view was 2°. 6 x 2°. 6 
and in the OCDM cosmology, 2°. 8 x 2°. 8. The angular reso- 
lution in the LCDM cosmology was 0'.34, which equates to 
the minimum value of the particle softening at the optimum 
redshift, z = 0.36, for lensing of sources at a redshift of 1. 
In the case of the OCDM cosmology, the angular resolution 
was 0'.37. In both cosmologies, to allow for the larger angu- 
lar size of the minimum softening at low redshifts and also 
for the range of particle softening scales above the minimum 
value, we have not declared our results below angular scales 
of 1 arcminute. 

To obtain the lensing statistics for sources at redshift 
z s — 1, we computed the required values assuming the 
sources were located at the front face of the simulation vol- 
ume whose redshift was closest to 1. Therefore in the LCDM 
cosmology, the sources were at z s = 0.99 and in the OCDM 
cosmology they were at z s = 1.03. 

Once the 3-d shear had been computed at each location, 
the multiple lens-plane theory, described clearly in Schnei- 
der, Ehlers & Falco (1992), was applied along the lines of 
sight. The procedure involved combining the intermediate 
Jacobian matrices computed using the appropriate values 
of the angular diameter distances in the relevant cosmol- 
ogy. Finally, the Jacobian matrices appropriate at z = for 
sources at z s = 1 for every line of sight were obtained. 

From the Jacobian matrices in each of the 10 simula- 
tion runs for each cosmology, we computed the shear com- 
ponents, 71 and 72, and the full shear, 7. The values for the 
smoothed components were computed on angular scales of 
l'.O, 2'.0, 4'.0 and 8'.0, using a top-hat filter. 

The computed values for the probability distributions 
from each of the N = 10 runs in each cosmology were aver- 
aged, and the errors on the means of la/VN determined. 



6 COMPARISON OF ANALYTICAL 

PREDICTIONS WITH RESULTS FROM 
NUMERICAL SIMULATIONS 

We compare here our analytical predictions for the smoothed 
shear components and for the shear modulus with the results 
of the numerical simulations. The quantites plotted are for 
various smoothing angles from 1' to 8' in the two different 
cosmologies, LCDM and OCDM, and in each case the source 
redshift is fixed at unity. 

Firstly we compare the results for the shear compo- 
nents, plotted in Figure 2, for the LCDM cosmology, and 
Figure 3, for the OCDM cosmology. Since the PDFs for the 
components are even functions we only plot the PDF, "P(7i), 
over 71 > 0. The solid lines correspond to our analytical 
predictions, (62) and (64), based on the stellar- model. For 
comparison, we also plot the Gaussian distribution which 
has the same variance (the dashed curves), so as to clearly 



LCDM 
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Figure 6. The PDFs, V(~fi) and V(\-f\), on a linear-scale for 
the LCDM cosmology and smoothing angle 9 S = 1'. The source 
redshift is again fixed at unity. The curves and the data-points 
are the same as in the upper left panels of Figures 2 and 4. 



show the deviations from Gaussianity (this corresponds to 
a Gaussian underlying density field). We find that at small 
angular scales, 9 S < 4', the results from the numerical simu- 
lations are very accurately described by our non-linear pre- 
dictions where there are clear signs of departure from Gaus- 
sianity, especially in the tails of the PDF. Of course, this 
merely reflects the non-linearity of the underlying mass dis- 
tribution probed by these small angular windows. At larger 
smoothing angles, 8 S > 4', the stellar Ansatz seems to over- 
predict the non-linearity induced by gravitational clustering 
and the results lie between the fully non-linear predictions 
and the Gaussian approximations. This behaviour may be 
due to the simple interpolation of eq.(47) for the transi- 
tion between the linear and highly non-linear regimes, which 
may not be sufficiently accurate. However, we did not try in 
this work to improve over the simple interpolation (47) be- 
cause we expect the stellar model to be less accurate in the 
weakly non- linear regime. Indeed, in the quasi-linear regime 
the angular dependence of the many-body correlations can 
be derived from perturbation theory and one can see that for 
the three-point function it already violates eqs.(38)-(41). On 
the other hand, in the highly non-linear regime, the angu- 
lar dependence of the many-body correlations gets shallower 
and eqs.(38)-(41) actually agree with numerical simulations 
for the three-point function (Scoccimarro & Frieman 1999). 
Of course, at very large scales (which only probe the linear 
regime) our results would become exact (within the small 
angle approximation) as we would recover the simple Gaus- 
sian. 

The PDFs for the shear modulus are presented in Figure 
4, for the LCDM cosmology, and Figure 5, for the OCDM 
cosmology. The solid lines show our analytical predictions 
(68)-(74)-(75) from the stellar-model. For comparison, we 
also plot (dotted-dashed curves) the results obtained by as- 
suming that the two shear components are uncorrelated but 
their own PDFs are still given by our non-linear prediction; 
this corresponds to eqs.(76)-(77). Therefore, we see from Fig- 
ures 4 and 5 that the correlation between both components, 
71 and 72, only has a rather small influence on the PDF 
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of the modulus, V(\j\). For reference we have also shown 
the results for the Gaussian case (dashed curves), where 
the PDF for I7 is given by a Rayleigh distribution with 
a variance equal to that of the shear components (dashed 
curves). Comparison of our predictions with the results of 
the numerical simulations shows that there is a remarkable 
agreement in both cosmologies (as found for the shear com- 
ponents), especially at the smaller smoothing scales. At the 
larger scales, the agreement is less good for the same rea- 
sons as described above. The sample variance arising from 
the various realisations also makes the comparison more dif- 
ficult; as the size of the smoothing radius is increased, the 
number of patches with completely independent information 
decreases, thereby increasing the sample variance. In addi- 
tion, although we made use of 10 different simulation runs 
(as described in Section 5), they were all based upon the 
same iV-body realisation in each cosmology, so the separate 
runs may not contain completely independent information. 
Much larger simulations would be required to quantify any 
departure from the analytical predictions. 

Finally, we display in Figure 6 our results over a linear- 
scale for V rather than the logarithmic scale of the previous 
Figures, for the LCDM cosmology and the smoothing angle 
9 S = 1'. While the previous Figures emphasized the tails 
of the PDFs we can clearly see in Figure 6 the shapes of 
the PDFs near their maxima. In agreement with the dis- 
cussion above, we can check that our predictions match the 
deviations from Gaussianity observed in the numerical sim- 
ulations. Thus, the PDF, "P(7i), of the shear component, 
71, is more sharply peaked around 71 = than a Gaussian. 
Moreover, the location of the maximum of the PDF, P(|7|), 
of the shear modulus is accurately reproduced by our an- 
alytical model. These features might be useful, in addition 
to the extended tails seen in the previous Figures, to mea- 
sure in weak lensing surveys the departures from Gaussian- 
ity brought by the non-linear gravitational dynamics. This 
may also help to break the degeneracies in the measurement 
of the cosmological parameters. 



7 DISCUSSION 

Extending earlier studies we have developed analytical tech- 
niques based on a simple hierarchical Ansatz, which we 
called the stellar model, to compute various statistical quan- 
tities associated with cosmic shear. Thus, we have com- 
puted the full probability distribution functions associated 
with the smoothed shear components, V{^a{0 s )), and the 
modulus of cosmic shear, T'(\~f(9 s )\). Here |7(# s )| is defined 
by first smoothing both shear components and second tak- 
ing their modulus. We have also derived the joint-PDF, 
r P(~/i(9s),"f2{6s))- These results extend similar studies ap- 
plied to the convergence field. However, being a spin-2 field, 
statistics from the shear are much richer than those from 
the convergence field and, more importantly, can be di- 
rectly used for observational purposes where the non-trivial 
survey geometry can pose a considerable problem for re- 
constructing convergence n(9 s ) maps. 

The comparison of our analytical calculations with the 
results of the numerical simulations shows that our simple 
method yields accurate predictions at small angular scales, 
9 S < 4' (for a source redshift z s — 1) which probe the non- 



linear regime for the underlying density field. Therefore, we 
can use with confidence the simple stellar model we intro- 
duced in this paper to study the shear statistics and ex- 
tract useful information from future surveys. In particular, 
our analysis shows how the deviations from Gaussianity in 
the density field built by the non-linear gravitational dy- 
namics translate into the PDF of the shear. This can also 
be seen in the Figures where we drew Gaussian PDFs for 
comparison. Hence we may hope to derive from shear maps 
some key properties of the structure of the density field (e.g., 
the amplitude of the many-body correlations and the shape 
of the PDF of the smoothed density contrast) which could 
shed some light onto the non-linear regime of gravitational 
clustering. In this respect, note that the separation of the 
shear, 7, into the factor K m in and the "normalized shear," 
7, also provides a good separation between the dependence 
on the cosmological parameters and the projection effects 
(e.g., f2 m ,f2A and z s ) versus the amplitude and the non- 
Gaussianities of the fluctuations of the density field which 
characterize the collisionless gravitational dynamics itself. 
This is obvious for the convergence, k, as recalled at the end 
of Section 3.2, but this feature also holds for the shear, 7. 

However, in the near future, where the observational 
noise may still be important, one may first focus on the es- 
timation of the cosmological parameters. To this order, our 
study is complementary to present analysis of the first few 
higher-order correlation functions of the shear field (Zaldar- 
riaga & Scoccimarro, 2003, Schneider & Lombardi, 2003, 
Bernardeau, Mellier & van Waerbeke, 2002, and Takada & 
Jain 2003) based on results derived for comological spin-2 
fields such as CMB polarization, where shape dependence 
and scale dependence of the shear is used to study non- 
Gaussianity induced by the underlying mass distribution. In 
any case, a few measures of non-Gaussianity (e.g., the skew- 
ness of the convergence, k, see Pen et al., 2003) are useful 
to remove degeneracies between the various cosmological pa- 
rameters (e.g., Q m vs as) since a Gaussian is defined by only 
one quantity (its variance), while we would like to measure 
several quantities. Therefore, one may use our prediction for 
the PDF, Pi'-/), to build some useful new observables. 

Before we can use our method to extract some infor- 
mation from weak gravitational lensing surveys, there are 
some points which need to be addressed. In particular, we 
should take into account the effects of a finite width for the 
distribution of sources over redshift, as in actual surveys. It 
would also be desirable to include realistic noise, such as that 
arising from the intrinsic ellipticity distribution of galaxies, 
and Poisson sampling due to the discreteness of the source 
distribution. The finite size of the catalogue is also worthy 
of consideration. These effects are beyond the scope of this 
paper, where we have mainly investigated the accuracy of 
our approach to describe the physics of weak gravitational 
lensing by large-scale structures. Having shown here that 
our simple stellar model provides a useful basis, we plan to 
investigate these points in detail in future work. 

To incorporate noise in our calculations we need to con- 
volve the pdf of shear which we have derived here with the 
noise pdf which is generally assumed Gaussian. A complete 
analysis is beyond the scope of the present paper. However, 
an order of magnitude estimate of the signal to noise ra- 
tio for the kurtosis of the shear (lowest order moment) has 
been given by Takada & Jain (2002) . They conclude that for 
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reasonable values of survey parameters (rms intrinsic ellip- 
ticy cr e = 0.4, number density of source galaxies n ga i = 
30 arcmin~ and the survey area Qsurvey — 25 degree ) 
it is possible to achieve a high signal to noise ratio for 
the small angular scales which we have probed (although 
our cosmological models differ slightly with the model they 
have considered but such estimates are largely independent 
of the background cosmology). Since deviation from non- 
Gaussianity in the case of shear pdf is dominated by kurtosis 
one would expect roughly similar conclusions for the entire 
pdf too. However a detailed analysis in line with Munshi & 
Coles (2002) will be presented in a future work. 

On the other hand, the present study already gives some 
information about delicate points which are of theoretical 
and practical interest. In particular, the underlying approx- 
imation in most analytical calculations (including ours) is 
the Born approximation. It neglects higher order correction 
terms in the photon propagation equation. The error intro- 
duced by such an approximation in the quasi-linear regime 
has been studied by several authors (e.g. Schneider et al. 
1998). Although perturbative calculations tend to show that 
this error is negligible for lower order cumulants, clearly such 
an analysis is not possible in the highly non-linear regime. 
However, previous studies of convergence statistics and our 
present results regarding the shear PDF, which exhibit a 
good match between theoretical predictions and the simu- 
lation results, indicate that such corrections are negligible 
even in the highly non-linear regime. 

With regard to numerical simulations, previous work for 
the predicted convergence statistics were made in compar- 
ison with ray-tracing experiments. In the present case, our 
analytical results for the shear, which are complimentary 
to the convergence results, are compared successfully with 
the numerical computations based on computing the full 3- 
d shear along every line of sight. The detailed match of our 
various predictions not only verifies our analytical calcula- 
tions but also increases our confidence in the two different 
numerical methods. 
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